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The observed mass pattern of scalar resonances below 1 GeV gives preference to the tetraquark as- 
signment over the conventional qq assignment for these states. We present a search for tetraquarks 
with isospins 0, 1 /2, 1 in lattice QCD, where isospin channels 1 /2 and 1 have not been studied 
before. Our simulation uses Chirally Improved fermions on quenched gauge configurations. We 
determine three energy levels for each isospin using the variational method. The ground state 
is consistent with the scattering state, while the two excited states have energy above 2 GeV. 
Therefore we find no indication for light tetraquarks at our range of pion masses 344 — 576 MeV. 
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1. Introduction 

The observed mass pattern of scalar mesons below 1 GeV, illustrated in Fig. does not 
agree with the expectations for the conventional qq nonet. The observed ordering thk < "1^0(980) 
can not be reconciled with the conventional iis and Ud states since mus > mad is expected due to 
rus > m^. This is the key observation which points to the tetraquark interpretation, where light 
scalar tetraquark resonances may be formed by combining a "good" scalar diquark 

[qQ]a = £ahc[qlCY5Qc-QlCr5qc] (color and flavor anti - triplet) (1.1) 

with a "good" scalar anti-diquark [qQ\a The states [qq]i,^ 3^ i^'lhf.Sc form a flavor nonet of 
color-singlet scalar states, which are expected to be light. In this case, the 7 = 1 state [us] [ds] with 
additional valence pair ss is naturally heavier than the 7=1/2 state [M(jf][<i!s] and the resemblance 
with the observed spectrum speaks for itself. 

Light scalar tetraquarks have been extensively studied in phenomenological models [|]], but 
there have been only few lattice simulations [§, ^, 0]. The main obstacle for identifying possi- 
ble tetraquarks on the lattice is the presence of the scattering contributions in the correlators. All 
previous simulations considered only 7 = and a single correlator, which makes it difficult to dis- 
entangle tetraquarks from the scattering. The strongest claim for a as tetraquark was obtained for 
ntjc ~ 180 — 300 MeV by analyzing a single correlator using the sequential empirical Bayes method 
[^]. This result needs confirmation using a different method (for example the variational method 
used here) before one can claim the existence of light tetraquarks on the lattice with confidence. 

We study the whole flavor pattern with 7 = 0, 1/2, 1 and our goal is to find out whether there 
are any tetraquark states on the lattice, which could be identified with observed resonances a (600), 
k:(800) and a()(980). Our methodology and results are explained in more detail in 
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Figure 1: Observed spectrum of scalar mesons below 1 GeV (left), together with the expected spectrum for 
the nonet of scalar tetraquarks (middle) and compared with a typical qq spectrum (right). 



2. Lattice simulation 

In our simulation, tetraquarks are created and annihilated by diquark anti-diquark interpolators 

ff'=o = [ud][ad] , ^^='/^ = [ud][ds] , ^^=1 = [us][ds] . (2.1) 
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In each flavor channel we use three different shapes of interpolators at the source and the sink 

K = [qnQn] WM'n] , = WM'J , ^3 = M.] [qim ■ (2-2) 

Here qn and denote "narrow" and "wide" Jacobi-smeared quarks. We use exactly the same two 
smearings as applied in [|6|], which have approximately Gaussian shape and a width of a few lattice 
spacings. In order to extract energies E„ of the tetraquark system, we compute the 3x3 correlation 
matrix for each isospin 

=I^'^^'(0|^/(x,0^f (0,0)10)^^^3 = X(0|^/|«)(«|^f|0) ^-^"^ = Iw;/.-^"' . (2.3) 

X n n 

Like all previous tetraquark simulations, we use the quenched approximation and discard the 
disconnected quark contractions. These two approximations allow a definite quark assignment to 
the states and discard [qq] [qq] ^ qq ^ vac mixing, so there is even a good excuse to use them 
in these pioneering studies. We work on two^ volumes V = x T = 16^ x 32 and 12^' x 24 at 
the same lattice spacing a = 0.148 fm [Q]. The quark propagators are computed from the Chirally 
Improved Dirac operator [Q] with periodic boundary conditions in space and anti-periodic boundary 
conditions in time. We use mia = m„ d<3 = 0.02, 0.04 and 0.06 corresponding to = 344, 475 and 
576 MeV, respectively. The strange quark mass m^a = 0.08 is fixed from m^. The analysis requires 
the knowledge of the kaon masses, which are 528, 576, 620 MeV for mia = 0.02, 0.04, 0.06. 

The extraction of the energies from the correlation functions using a multi-exponential fit 
Cij = Y^n^n^^^"' is unstable. A powerful method to extract excited state energies is the variational 
method, so we determine the eigenvalues and eigenvectors from the hermitian 3x3 matrix C{t)^ 

C{t)v„{t)=Xn{t)Mt) ■ (2.4) 

The resulting large-time dependence of the eigenvalues A„(?) = w„e^^"' [1 + (^^{e^^"')] allows a 
determination of energies £0,1,2 and spectral weights wo,i.2- The eigenvectors v„(?) are orthogonal 
and represent the components of physical states in terms of variational basis ( |2.2[ ). 

3. Results 

Our interpolators couple to the tetraquarks, if these exist, but they also unavoidably couple to 
the scattering states nn (I = 0), Kn (/ = 1 /2) and KK, nrjss (t = 1) as well as to the heavier states 
with the same quantum numbers. The lowest few energy levels of the scattering states Pi (k)P2{—k) 

^^^^+^^2, ... , ^m2^+(^)' + ^m2,+ (^y,... (3.1) 

are well separated for our L and we have to identify them before attributing any energy levels 
E ~ m(,,K,ao to the tetraquarks. 

The three energy levels for 7 = are represented by the effective masses of the three eigen- 
values A„(?) in Fig. |2[ Similar effective masses are obtained for all isospins, quark masses and 

'Equal smearings are used on both volumes. 

^We use the standard eigenvalue problem in order to study w{L). 
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volumes. The energies £"0,1,2 extracted from Ao,i,2(0 are summarized in Fig. ^ for all isospin 
channels. The excited energies were obtained from a conventional two-parameter fit Ai,2(0 = 
>vi,2[^"^'^' + ^^''^^^^']- The fit of Ao(f) ( |A.1| ) takes into account a non-standard time dependence 
at finite T (noticed from decreasing nieff near f ~ r/2) and is described in Appendix. 
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Figure 2: Effective masses for the three eigenvalues Ao.i.2(f ) at / = 0, m/a — 0.04 and V — 16^ x 32. They 
are obtained from conventional cosh-type definition, except for the full diamonds which are obtained from 
Ao(f) - const = wo[e""''-/^' +e""V/(^-0] (see Appendix). The lines give energy levels for the scattering 
states: full lines present non-interacting energies, while dashed lines take into account tree-level energy 
shifts. 

The ground state energies in / = 0, 1/2 and 1 channels are close to Iniji, nij^+mK and Ithk, m^i + 
m,j^^ , respectively, which indicates that all ground states correspond to the scattering states Pi (0)^2(0) 
(see Fig. ^). Another indication in favor of this interpretation comes from wo(L = \2)/wq{L = 
16) ~ 16V 12^ (see Fig. |). This agrees with the expected dependence wq'x Xjl} for scattering 
states [Q], which follows from the integral over the loop momenta / j0jjf(k,t) — > jiJLlf(^:^) 
with dki = 2n/L. The third indication comes from a non-standard time dependence of Xo{t), which 
agrees with the expected time-dependence of Pi(0)P2{0) at finite T (see Appendix). 

The most important feature of the spectrum in Fig. ^ is a large gap above the ground state: 
the first and the second excited states appear only at energies above 2 GeV. Whatever the nature of 
these two excited states are, they are much too heavy to correspond to a (600), K"(800) or (20(980), 
which are the light tetraquark candidates we are after. The two excited states may correspond to 
P\ (k)P2{—k) with higher k or to some other energetic state. We refrain from identifying the excited 
states with certain physical objects as such massive states are not a focus of our present study. 

Our conclusion is that we find no evidence for light tetraquarks at our range of pion masses 
344 — 576 MeV. This is not in conflict with the simulation of the Kentucky group ||3|], which finds 
indication for an / = tetraquark for pion masses 1 80 — 300 MeV (but not above that) since all our 
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Figure 3: Three lowest energy levels from tetraquark correlators in / — 0, 1 /2, 1 channels at lattice vol- 
ume 16^ X 32. The lines give analytic energy levels for scattering states: full lines present non-interacting 
energies, while dashed lines take into account tree-level energy shifts. 



pion masses are just above 300 MeV. 

4. The absence of scattering states with \k\ = In /Lin the spectrum 

Why are there no states close to the energies of Pi(l)P2(— 1) with \k\ = In/L in the spectrum 
of Figs. ^ and ^? We believe this is due to the fact that all our sources CLl) have a small spatial 
extent and behave close to point like. The point source couples to all the scattering states equally 
Y,ig-i.\P\(Ji)P2i—k)) up to a factor gj, which gives a Lorentz structure of the coupling. In this 



approximation each source ^1^2,3 (2-2) couples to the few lowest scattering states equally (within 
our error bars) 

ffi = CiY,gl\Pi(k)P2i-k))+ai\a)+bi\b) / = 1,2,3 (4.1) 

1 

and only the overall strengths ci 2,3 are different. We assume also that the sources couple to two 
other physical states a and b. Given these linear combinations, one can construct the corresponding 
3x3 correlation matrix and it can be easily shown that its three eigenvalues are 



,-(£nft+£'=2(J))r 
EPi{k)EP2{k) 



■ Wa.b e 



-Ea,bt 



(4.2) 



The physical states a and b get their own exponentially falling eigenvalues, while a tower of few 
lowest scattering states contributes to a single eigenvalue in this approximation. The agreement 
of Ao(f) and the analytic prediction {X2) with = 1 [§1 gives us confidence about the hypothesis 



5 



A lattice study of light scalar tetraquarks with isopins 0, 1/2 and 1 



Sasa Prelovsek 



C 

o 

.o 



1=0 



3 3 

16/12 



0.1 



0.2 



0.3 



m [GeV] 



1=1/2 



0.1 



I I \ 

0.2 0.3 



Figure 4: The ratio of spectral weights wo(L = \2)/wq{L — 16) for / = 0, 1/2 as computed from ground 
state eigenvalues for two volumes Lr' . 

So our basis does not dissentangle the few lowest scattering states into seperate eigenvalues. 
They would contribute to seperate eigenvalues only by using a different or larger basis. 

5. Conclusions and outlook 

Our lattice simulation gives no indication that the observed resonances a (600), k:(800) and 
ao(980) are tetraquarks. However, one should not give up hopes for finding these interesting objects 
on the lattice. Indeed, our simulation with pion masses 344 — 576 MeV does not exclude the 
possibility of finding tetraquarks for lighter niji or for a different interpolator basis. A stimulating 
lattice indication for a as a tetraquark at m-n = 182 — 300 MeV has already been presented in [^. 
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A. Appendix: Effect of finite T on scattering states 

We find that the cosh-type effective mass for the ground state Pi (0)^2(0) is decreasing near t ~ 
r/2 (empty symbols in Figs. ^ and Q), which means that Ao(f) does not behave as g^^of _|_g-£o(r-0 
at large t. The time dependence of Pi(0)P2(0) with anti-periodic propagators in time is 

^Pi(o)P2(o)^^^ = WQ[e-E^' +e-'^'^^-'^] + A +e-'««fg-m/.i(T-0] ^ ~ mpi +mp2 . 

(A.l) 
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In the last term one pseudoscalar propagates forward and the other backward in time (see Appendix 
A of [^). The ground state energies for isospins and 1/2 in Fig. ^ were obtained by fitting 
to (AJ_) with three unknown parameters £'o> and A ^ . The effective mass obtained from 



h){t) -A[£'-'«^iff'-'"«(^-') (^-')] is flat near f ~ r/2 (full symbols in Figs. | andl). 

The 7 = 1 channel contains KK and tttj^ scattering states and is therefore more challenging. 
The ground state energies were obtained from a naive one-state (cosh) fit since a two-state fit is 
very unstable. 
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Figure 5: The effective masses for the 1—1/2 ground states at various m/o and V = 16^ x 32. Empty 
symbols are obtained using the conventional cosh-type effective mass Ao(r) — vvo[e^"''''^' + e^'"<-//(^^')], 
while the full symbols take into account the correct form Ao(f) - A[e""''''e""''^(^"') ^ g-mKtg-m^{T-r^ ^ 
wole^'^^'/f' +e'"''ff'-'^^'^] for/:(0)7r(0) scattering (see Appendix). 
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At the time of our presentation at Lattice 2008, we were not aware of the last term in ( A. 1), so the numerical results 



on our slides are slightly different. The general physical conclusions concerning the tetraquarks are however the same. 
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